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ABSTRACT 

Aims. We determine the components of the A-effect tensor that quantifies the contributions to the turbulent momentum 
transport even for uniform rotation. 

Methods. Three-dimensional numerical simulations are used to study turbulent transport in triply periodic cubes under 
the influence of rotation and anisotropic forcing. Comparison is made with analytical results obtained via the so-called 
minimal tau-approximation. 

Results. In the case where the turbulence intensity in the vertical direction dominates, the vertical stress is always 
negative. This situation is expected to occur in stellar convection zones. The horizontal component of the stress is weaker 
and exhibits a maximum at latitude 30° — regardless of how rapid the rotation is. The minimal tau-approximation 
captures many of the qualitative features of the numerical results, provided the relaxation time tau is close to the 
turnover time, i.e. the Strouhal number is of order unity. 
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1. Introduction 

Differential rotation plays a crucial role in dynamo pro- 
cesses that sustain large-scale magnetic activity in stars like 
the Sun (e.g. Moffatt [Mi Krause & RadlerpJHOl)- The in- 
ternal rotation of the Sun is familiar from helioseismology 
(e.g. Thompson et al. 2003|), but the processes sustaining 
the observed rotation profile are not understood well. The 
angular momentum balance of a rotating star is determined 
by the conservation equation 
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where U is the meridional flow, s the cylindrical radius, p 
the density (neglecting however its fluctuations), Q = U^/s 
the local angular velocity, and u$u the zonal component 
of the Reynolds stress. Overbars denote averages over the 
azimuthal direction. 

The meridional flow can also be directly affected by 
the Reynolds stresses (e.g. Riidiger Il989p . but it is more 
strongly determined by the baroclinic term that arises if the 
isocontours of density and pressure do not coincide. Such 
a configuration can appear because of latitude-dependent 
turbulent heat fluxes that arise naturally in rotating con- 
vection (e.g. Riidiger [1982; Pulkkinen et al. 1993; Kapyla 
et al. 120041 Riidiger et al. I2005ap or from a subadiabatic 
tachocline (Rempel I20U5|) which is likely to occur in the 
Sun (Rempel HOni Kapyla et al. I2006bjl . The flows due to 
thermodynamic effects are likely to be needed to avoid the 
Taylor-Proudman balance in the Sun (e.g. Durney 1989: 
Brande nburg et al. [T9"9"2l Kitchatinov & Riidiger fUMA 
Rempel 2005). The overall importance of the meridional 



flow in the angular momentum balance of the Sun is, how- 
ever, still unclear since no definite observational informa- 
tion about it is available below roughly 20 Mm depth (e.g. 
Zhao k, Kosovichev l2"UM[) . 

Although not much more is known about the Reynolds 
stresses from observations, already this limited knowledge 
can be used to gain insight into the theory of turbulent 
transport. Solar surface observations indicate that there is 
an equatorward flux of angular momentum, as described 
by the Reynolds stress component Qg^ = uqu^, of sev- 
eral 10 3 m 2 s~ 2 in the latitude range where sunspots are 
observable (e.g. Ward 1965; Pulkkinen & Tuominen 1 1998; 
Stix[2002). In mean-field theory the simplest approximation 
that can be made concerning the Reynolds stresses is to as- 
sume them proportional to the gradient of mean velocity 
(the Boussinesq ansatz): 

Qij = muj = -AfijkiUkj ■ (2) 

In the Sun this ansatz turns out to be insufficient because 
the observed Qg^, and solar surface differential rotation pro- 
file indicate that the turbulent viscosity is negative. Thus, 
in analogy to mean-field dynamo theory, additional contri- 
butions to the Reynolds stress were conjectured to appear 
(e.g. Wasiutynski 1946; Krause & Riidiger II 974ft . leading 
to the present description 

Qij = Aijkflk — NijkiU k,i , (3) 

where Kijk is a third-rank tensor describing the A-effect 
that contributes to the Reynolds stress even in the case of 
rigid rotation. These terms are often referred to as "non- 
diffusive" contributions to the Reynolds stress. The zonal 
components of the stress can be written in the form (e.g. 
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Ay sin 8 Q — (1 — e)r sin 9 



dVL 

dr 



(5) 



where Ah and Ay describe the non-diffusive transport and 
where v t is the turbulent viscosity. The factor 1 — e in the 
latter equation indicates that the turbulent viscosity can 
be anisotropic. Furthermore, the coefficients Ah, Ay, and 
v t can vary as functions of the spatial coordinates. 

Much effort has been devoted to determining Reynolds 
stresses from convection simulations (Hathaway & 
Somerville 119831 Pulkkinen et al. 1993; Rieutord et al. 
[19941 Brummel l et al. [1993 Chan [20 011 Kapyla e t al. [20041 
Rudi ger et al. I2005bl Hupfer et al. I27KI51 l27K?r31 Giesecke 
2007). These studies have confirmed the existence of the 
A-cffect and revealed some surprising results that are at 
odds with theoretical considerations (e.g. Kitchatinov & 
Rudiger 1993, 12005P derived under the second-order corre- 
lation approximation (SOCA). The discrepancies are most 
prominent in the rapid rotation regime, f2* ~ 10, where 



= 2 f^oTto 



(G) 



is the Coriolis (or the inverse Rossby) number. Here, f2o 
is the angular momentum-averaged rotation rate and r to 
the convective turnover time. In the solar convection zone, 

varies between 10 -3 near the surface to ten or more in 
the deep layers. Convection simulations in the latter regime 
show that the horizontal angular momentum flux is directed 
toward the equator, corresponding to positive Qg^ in the 
northern hemisphere, and that it peaks very sharply near 
the e quator (Chan [10011 Kapyla et al. [HI Hupfer et al. 
2005). On the other hand, the vertical stress Q r( f, can be 
directed outward (Kapyla et al. I2004p . contradicting the 
theory for vertically dominated turbulence (e.g. Biermann 
[TMT1 Riidiger fTUSUl [1^591) . So far, these results remain with- 
out proper explanation. 

Often the Reynolds stress realized in the simulation is 
taken to solely represent the A-effect. This approach seems 
like a reasonable starting point but in an inhomogeneous 
system large scale mean flows are generated when rota- 
tion becomes important. These flows affect the Reynolds 
stresses via the turbulent viscosity. Furthermore, in the 
presence of stratification, heat fluxes can also significantly 
affect the stresses (Kleeorin & Rogachevskii 200|)|). In the 
present study we simplify the situation as much as possible 
in order to disentangle the effect of the turbulent veloc- 
ity field from other effects. Thus we neglect stratification 
and heat fluxes by adopting a periodic isothermal setup. 
Turbulence is driven by external forcing, which provides 
clear scale separation between the turbulent eddies and the 
system size, which is typically not achieved in convection 
simulations. Further insight is sought from comparison of 
simple analytical closure models with numerical data. 

Preliminary results on the A-effect are presented in 
Kapyla & Brandenburg (2007). In the present paper, nu- 
merical datasets covering a significantly larger part of the 
parameter space are analyzed, and a more thorough study 
of the validity and results of the minimal tau-approximation 
are presented. 

The remainder of the paper is organized as follows. 
Section [U summarizes the model and the methods of the 
study, and in Sects. [3] and [4] the results and the conclusions 
are given. 



2. The model and methods 

2.1. Basic equations 

We model compressible hydrodynamic turbulence in a 
triply periodic cube of size (27r) 3 . The gas obeys an isother- 
mal equation of state characterized by a constant speed of 
sound, c s . The continuity and Navier-Stokes equations read 
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(8) 



where D/Dt = d/dt + U ■ V denotes the advective deriva- 
tive, U is the velocity, p the density, / v isc the viscous force, 
and /force the forcing function. Due to the periodic bound- 
aries, mass is conserved and the average density retains its 
initial value ~p = po at all times. Compressibility is retained 
but we consider low Mach number flows, w rms /c s « 0.1—0.2. 
The viscous force is given by 



/vise = ^(V 2 £/+±VV-[/ + 2S- Vlnp 



where v is the kinematic viscosity and 



_ 1 (dU i dUj , 1 
oxj oa 



dU k 
* UlJ dx k 



(9) 



(10) 



the traceless rate of strain tensor. The forcing function 
/force is given by 



f{x, t) = Re{AT • f k(t) exp[ifc(t) • x - i<p(t)}} 



(11) 



where x is the position vector, TV = /(°) 'c s (fcc s / 'St) 1 / 2 is 
a normalization factor, and describes the direction- 
dependent amplitude of the forcing (see below), k = \k\, St 
is the length of the time step, and — 7r < <j){t) < ir a random 
delta-correlated phase. The vector /& is given by 



fk = 



k x e 



y/k 2 - (fe • e) 2 



(12) 



where e is an arbitrary unit vector. Thus, fk describes non- 
helical transversal waves with \fk\ 2 = 1 , where k is chosen 
randomly from a predefined range in the vicinity of the av- 
erage wavenumber kf/ki = 5 at each time step, where fci is 
the wavenumber corresponding to the domain size, and fcf 
the wavenumber of the energy-carrying scale. 

To make the resulting turbulence anisotropic the forcing 
amplitude depends on direction 



(<%/o 



cos 2 e fc /i)/ 



(iso) 



(13) 



where fo and /i are the amplitudes of the isotropic and 
anisotropic parts of the forcing, z the unit vector in the 
vertical direction, and 0& the angle between the vertical 
direction and the wave vector k. In the following we as- 
sume fi 3> fo, but even then the amount of anisotropy 
depends on the Reynolds number and is in any case only 
quite modest (see, e.g. Table iBTj) . 

At this point it is important to emphasize that our forc- 
ing function is designed to capture the effects that lead to 
a finite A-effect. The implementation of anisotropy should 
therefore be a simple one. In stars, anisotropy is produced 
by stratification and convection. Our goal is clearly not to 
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simulate properties of convection other than its tendency 
to make the turbulence anisotropic. 

The numerical computations were made with the 
Pencil-CodeQ, which uses sixth-order accurate finite dif- 
ferences in space, and a third-order accurate time-stepping 
scheme (Brandenburg & Dobler [20021 Brandenburg 2003). 
Resolutions up to 256 3 grid points were used in the simu- 
lations. 



2.2. Nondimensional units 

In the following we use non-dimensional variables by setting 

c„ = fci = po = 1- (14) 
This means that the units of length, time, and density arc 



[x]=ki\ [t] = (c s ki)- 1 , [p]=p . 



(15) 



However, in most of the plots we present the results in 
explicitly non-dimensional form using the quantities above. 

2.3. Coordinate system, averaging, and error estimates 

The simulated domain is thought to represent a small rect- 
angular portion of a spherical body of gas. We choose 
(x, y, z) to correspond to (9,<f>,r) of spherical coordinates. 
With this choice the rotation vector can be written as 



ft = f7 (-sin6»,0,cos6l) rj 



where 9 is the angle between the rotation axis and the local 
vertical direction, i.e. the colatitude. 

Since the turbulence is homogeneous, volume averages 
are employed and denoted by overbars. An additional time 
average over the statistically saturated state of the calcu- 
lation is also taken. We define the Coriolis number as 



Co = 



(17) 



where u r ms is the rms-velocity. In comparison to the com- 
monly used definition definition (TTTJ) is smaller by a fac- 
tor of 27r, i.e. Co = i\/27r. The other dimensionless number 
relevant in this study is the Reynolds number based on the 
forcing scale 



Re = 



vkf 



(18) 



See Fig. Q] for a snapshot from a typical run with Re w 60. 

Errors are estimated by dividing the time series into 
three equally long parts and computing mean values for 
each part individually. The largest departure from the mean 
value computed for the whole time series is taken to repre- 
sent the error. 



2.4. The K-ef Feet from the minimal tau-approximation 

To have some understanding of the numerical results, 
we compare with the simplistic tau-approximation (here- 
after MTA) in real space (e.g. Blackman & Field |2"0D21 
Brandenburg et al. [5004). Unlike the usual first-order 
smoothing approximation where nonlincaritics in the fluc- 
tuations are neglected, they are retained in an approximate 
manner in MTA. 



http : //www.nordita. org/sof tware/pencil-code/ 



t= 200.0 
u 




Fig. 1. U x at the periphery of the simulation domain from 
a slowly rotating run with Co ~ 0.3, fcf/fci = 5, and = 0°, 
resolution 256 3 . 



In order to develop a theory for the Reynolds stress, 
Qij = UiUj, we derive an equation for its time derivative, 



(!6) = iiiUj + 



UiUj 



(19) 



In the absence of large-scale flows, i.e. U = 0, the equation 
for the fluctuating part can be written as 



(20) 

where Ni = —u^dkUi — c^c^hi/? is a nonlinear term. 
Multiplying Eq. (|20|) by Uj gives 



UjXli UjNi 2 Zlj£i 7nn fl rn U n ~\~ Ujfi 

Inserting this into Eq. (TT9"]) yields 



(21) 



-2 Cjki^lkQu — 2 eikiflkQji + u%fj + Ujfi + Tij , (22) 



where = UiNj +UjNi are the triple correlations. Under 
the assumption of periodic boundary conditions, this term 
can be written in the form 



= — c s (uidj In p + Ujdi In p) + UiUj V • u . 



(23) 



In MTA the higher than second-order terms, i.e. Tij, in 
the equations of turbulent correlations are retained in a 
collective manner by parametrizing them with a term that 
is equal to the original correlation divided by a relaxation 
time, i.e. 



Ti 



(24) 



Now the equations can be solved for the stresses in terms of 
the forcing. A simple ansatz for parameterizing the forcing 
is given in terms of the non-rotating equilibrium solution, 



Uifj +u J f l 



(25) 



Inserting the parameterizations ([2~4")l and (|25p into Eq. 
yields 



Qi 



-2 tjki^kQu — 2 €iki^kQ 



J I 



(0)s 



(26) 
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Fig. 2. Volume-averaged Reynolds stress components Q xx 
(top), Q yy (middle), and Q zz (bottom), normalized by the 
square of the rms-velocity, as functions of latitude and ro- 
tation from the turbulence simulations listed in Table IB. II 
Coriolis number, as defined in Eq. (| 1 T[) . varies as indicated 
in the legend in the middle panel. 



see Appendix for more details on the equations used in 
the MTA-model. The model equations are very similar to 
those of Ogilvie (|2003l see also Garaud & Ogilvie [2~005p who 
studied hydrodynamic and magnetohydrodynamic turbu- 
lence and angular momentum transport due to shear flows 
and the magnetorotational instability. Among other things, 
they also introduced an isotropization term that causes de- 
caying turbulence to become isotropic. Some of our models 
give explicit evidence of such a term. 



Fig. 3. Diagonal Reynolds stress components from the 
MTA-model. Here the Coriolis number is defined via 
Eq. ([23) with St = 1 and f = 0. 

The rotational influence in the MTA-model is measured 
by the Coriolis number 

CoMTA^^f^CoSt, (27) 

where 

St = ru lms k{ (28) 

is the Strouhal number, which is the main free parameter 
in the model. 

On account of previous results on similar systems of 
forced turbulence (e.g. Brandenburg & Subramanian 2005) 
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Fig. 4. Same as Fig. [3] but with St = £ = 1. Compare with 
Figs. [3 and G2 



St = 1 is used as our reference model. To reproduce the 
numerical results, we introduce an empirical isotropization 
term in the MTA-model; see Sect. 13.11 and in particular 
Eqs. (|29|) and (j30|) . The use of this term is regulated by 
another free parameter £, which can obtain the values zero 
or unity. 

The MTA-model results were obtained by advancing the 
time-dependent equations (|22p with the parameterizations 
(|2"4"|) and (|2"5)) until a stationary solution was reached. 



3. Results 

3.1. Diagonal components of the stress 

Consider turbulence where the intensity in the vertical di- 
rection is stronger than the intensities in the horizontal 
directions. This situation is encountered in the turbulence 
caused by the convective instability. First we consider the 
case with the maximum anisotropy that was achieved with 
the present model. Table IB. II summarizes the results for 
seven values of the Coriolis numbers ranging from 0.06 to 
roughly 5.4. Calculations at seven latitudes from the pole 
(9 = 0°) to the equator (9 = 90°) with equidistant inter- 
vals of A8 — 15° were made with each Coriolis number, 
Re « 12 ... 14 in all runs. In this set of runs, fo = 10~ 6 and 
fx = 0.2 were used in Eq. (TT3j) . 

Figure[5]shows the diagonal components of the Reynolds 
tensor as functions of latitude and Coriolis number from the 
numerical turbulence simulations. As rotation is increased, 
the magnitudes of the horizontal components Q xx and Qyy 
increase monotonically while Q zz decreases, which illus- 
trates an isotropizing effect of rotation on the turbulence. 
In the following, we refer to this effect as the rotational 
isotropization of turbulence. This effect is a purely empir- 
ical and refers to the observation that increased rotation 
leads to stronger mixing which washes out anisotropics that 
are caused by other effects. Of course, rotation itself can 
cause the turbulence to become anisotropic, but this seems 
to play a role only at much higher rotation rates. This is 
indeed seen in the most rapidly rotating case where the be- 
havior is more complex. The MTA-model, Fig. [31 on the 
other hand, fails to reproduce isotropization of turbulence 
when rotation is included. The behavior is most obvious at 
the pole where rotation contributes no net effect to linear 
order, see Appendix [Al Eqs. (|A.1[) to (|A.6|) . The same is 
true for the Q xx component at the equator. 

Some persistent trends arise as a function of latitude in 
the numerical simulations; for instance, Q xx peaks at the 
pole and decreases toward the equator except again for the 
fastest rotation. An approximately opposite trend is seen 
for Qyy for slow and rapid rotation, whereas in the interme- 
diate range a minimum appears at mid-latitudes. For slow 
rotation, Q zz behaves approximately in the opposite way to 
Qyy, having a maximum at the pole and a minimum at the 
equator, whereas for rapid rotation the trend is reversed. 
There is also a persistent maximum at mid-latitudes. The 
MTA-model, however, fails to reproduce most of the charac- 
teristics of the latitude distribution of the diagonal stresses. 
This indicates that nonlinear terms, which are manifestly 
not described adequately well by the MTA-assumption, are 
the deciding factor in determining the behavior of the di- 
agonal stresses. 

To capture the rotational isotropization with the MTA- 
model at least qualitatively, we experimented by adding a 
term 



Qi 



-^rot {Qij q $ij Q ) i 



(29) 



on the rhs of Eq. ([26]) . Here Q is the trace of Qij and 



3 Co 



MTA 



1 + Co 



(30) 



MTA 



The functional form of F rot is chosen purely empirically so 
that the magnitudes of the off-diagonal stresses in compari- 
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son to the simulations are fairly accurately reproduced with 
the MTA-model (see Fig. EJ. 

Figure 2] shows the results for the diagonal components 
with St = £ = 1. Now the magnitudes of the turbulence 
intensities are more in line with the full numerical sim- 
ulations, although the latitude distribution is still man- 
ifestly wrong. Although the off-diagonal components are 
much better represented by the linear terms appearing in 
the equation of the Reynolds stress (see Sect. I3.3|) . the ro- 
tational isotropization term helps for reducing their magni- 
tudes closer to the levels seen in the direct simulations also 
in that case. 

The functional form of F mt in Eq. (f2"9")) indicates that 
-Plot - > 3 for rapid rotation. This behaviour cannot be jus- 
tified based on the present numerical data (see Fig. [2]). 

3.2. Anisotropy of the turbulence 

Turbulence anisotropies can be characterized by the quan- 
tities 

A H = Qvy 7 Qxx , (31) 
it 

rms 

a _ Qxx + Qyy — 2 Q zz , . 

11. 

rms 

The importance of these quantities is that, for slow rota- 
tion, they can be considered as proxies of the A-effect ac- 
cording to Ah ~ 2t c Ah and Ay ~ 2t c Av (e.g. Riidiger 
1980, 1989), where r c is the correlation time of the turbu- 
lence. 

Table IB. II shows that for slow rotation Ah increases 
monotonically from the pole to the equator. For Co > 0.3, 
Ah exhibits a negative minimum at mid-latitudes and 
reaches a positive maximum at the equator. The maxi- 
mum at the equator can be explained as resulting from 
the surviving (y-)component of the Coriolis force at low 
latitudes. The relation between the horizontal A-effect and 
the anisotropy parameter Ah is at best poorly confirmed 
by the numerical results. It must also be noted that Ah is 
smaller by at least an order of magnitude in comparison 
to Ay, which also enters the equation for Q xy but with a 
higher order in Co (sec Eq. IA.7|) . 

The vertical anisotropy, Ay, on the other hand, re- 
tains its sign in all models. For slow rotation the absolute 
value of Ay decreases monotonically from the pole toward 
the equator. The latitude distribution for rapid rotation is 
approximately opposite, although an additional minimum 
(Co w 0.6 ... 1.5) or a maximum (Co ~ 5.4) can occur at 
mid-latitudes. Here the correspondence between Ay and 
Ay holds at least for the sign for all calculations, bar the 
few cases in the intermediate rotation regime where Q yz is 
positive (see below). 

3.3. Off-diagonal components of the stress 

In stellar convection zones, the off-diagonal Reynolds 
stresses contribute to the angular momentum transport and 
work to generate (A-effect) or to smooth out (turbulent vis- 
cosity) differential rotation. In the present case only the for- 
mer effect is in operation. Figure [3J summarizes the results 
for the runs listed in Table [B~T1 

The Q xy component of the stress corresponds to Qg^ in 
spherical coordinates and is responsible for latitudinal an- 



gular momentum transport. In the simulations this compo- 
nent is always positive and the latitude distribution peaks 
at latitude 30° (see the uppermost panel of Fig. [5]). The 
sign is in accordance with solar observations (Ward 1965; 
Pulkkinen & Tuominen [T998) and analytical turbulence 
models (Kitchatinov & Riidiger IT9931 12005)) . The latitude 
distribution in the rapid rotation regime is significantly dif- 
ferent from convection simulations where Q xy is sharply 
concentrated near the equator (Chan [200 1; Kapyla et al. 
120041 Hupfer et al. 120051 Riidiger et al. l2005ajl . The reason 
for this difference is still unclear but is evidently related to 
the physics that have been omitted in the present study. 
The MTA-model captures the latitude dependence rather 
well, but the magnitude of the stress is clearly too great 
(see Fig. [6]) . If the rotational isotropization term, Eq. (|29]) , 
is taken into account (see the uppermost panel of Fig. [7J , 
the agreement is also better for the magnitude. Keeping the 
forcing and rotation rate Qq fixed, the best agreement with 
the 3D models is found if St = 2 is used, see Fig. [8] 

While the analogue of the Q xz component does not play 
a direct role in the angular momentum balance in stars, it 
can still contribute via generating meridional flows (e.g. 
Riidiger 1989). In the numerical simulations we find that, 
for all cases except the most rapidly rotating one, Q xz is 
negative and peaks at 8 = 45° (middle panel of Fig. [5]). 
For Co w 5.4, however, the sign changes near the equator, 
with positive values toward the equator and negative ones 
toward the pole. This is at odds with the analytical result 
of Riidiger et al. ( 2005b]), but does agree with the results 
of convection simulations (Pulkkinen et al. I1993| Kapyla et 
al. 1200411 . 

The MTA-model gives qualitatively similar results, al- 
though the sign change occurs at significantly slower rota- 
tion; see the middle panel of Fig. [6] Rotational isotropiza- 
tion helps to correct the magnitude, but not the earlier 
occurrence of the sign change (Figs. [7] and [5]). 

Since the turbulent intensity of the vertical (z-)motions 
is greater than the horizontal ones, the expectation is that 
the stress component Q yz is negative, i.e. that Ay oc Ay 
(Biermann I195ip . This is indeed seen in the simulations 
quite consistently (lowermost panel of Fig. [5]), although at 
intermediate rotation low positive values can occur at high 
latitudes. The highest values of Q yz occur for Co i=s 0.3 
as opposed to Co w 0.6 for Q xy . Rotational quenching of 
Q yz seems to be stronger and occur for lower Co than for 
Qxy (see also Fig. (T0J. A similar trend was seen in convec- 
tion simulations by Kapyla et al. (2004). The consistently 
positive values of Q yz for rapid rotation seen in convection 
simulations (Kapyla et al. 20QU Chan 2007, private com- 
munication) do not occur in the present calculations. 

The value of Q yz always reaches a maximum at the 
equator. This contradicts with analytical results derived 
under first-order smoothing (Kitchatinov & Riidiger 1993, 
2005) and the MTA-model that predicts a maximum 
around 8 = 45° for intermediate and rapid rotation. The 
rotational isotropization term is again needed to reduce 
the magnitude of the stress. It seems that to reproduce 
Qyz correctly, one would need to apply somewhat stronger 
isotropization than what is presently used (see Figs. [7] and 

El- 

One conclusion that can be drawn from the MTA-model 
results is that, although the diagonal Reynolds stresses are 
quite poorly reproduced in comparison to the 3D simula- 
tions, the off-diagonals have most of the qualitative fea- 
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tures correct. This seems to imply that, to model the off- 
diagonals, the exact parameterization of the nonlinearities 
in the equation of the Reynolds stress is not crucial. The 
empirical rotational isotropization term helps to capture 
some of the missing features for the diagonal components 
and reduces the magnitudes of the off-diagonals to the level 
that is also seen in the simulations. 

3.3.1. Comparison with inhomogeneous simulations 

The homogeneous setup used so far prevents any mean 
flows from being generated. This is good for the purpose 



of testing the sole effect of turbulent velocity field on the 
Reynolds stresses, but can be argued to be unphysical be- 
cause astronomical objects where turbulence is important; 
e.g., the solar convection zone has boundaries and cannot 
be considered homogeneous. 

To test how much the assumption of homogeneity af- 
fects the results, we made a set of simulations with a setup 
where the z-boundaries are impenetrable. For the horizon- 
tal velocity components, stress- free boundary conditions 
are used, i.e. 

U x , z = U VtZ = U 2 = Q. (33) 
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Fig. 9. Same as Fig. [S] but for homogeneous (left panel) and inhomogeneous (right) simulations. 



To compare with the homogeneous simulations and at the The exception to this trend is the equator, i.e. — 90°, 
same time minimize the effects of the boundaries, we av- where a large mean shear flow develops, for which 



erage the results over — jir < z < j7T and the full (x,y)- fcjj 



extents. 



-fin 



(34) 



Similar flows have been seen in convection simulations 
(Chan HUE! Kapyla et al. 130031 Brandenburg [3007D . The 
volume average of this flow is zero. For the present simu- 
Mean flows are generated in the runs that, however, are lations the mean velocities near the z-boundaries can ap- 
small in comparison to the rms-velocity of the turbulence, proach the speed of sound or even exceed it. This prevents 
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us from computing models with f^o 
tistically saturated state. 



0.5 ... 2 until a sta- 



The results for the off-diagonal stresses are compared 
in Fig. [9] For slow rotation, Co < 0.62, the differences are 
minute at all other latitudes except at the equator. The 
tendency for Q yz to have a minimum at the equator is rem- 
iniscent of results from convection simulations (e.g. Chan 



I200T1 Kapyla et al. [20041 Rudiger et al. 12005b)) . For more 
rapid rotation, the trends for different components seem to 
diverge: Q xy is somewhat reduced whereas Q xz seems to 
increase somewhat and there is hardly any change for Q yz 
apart from the equatorial case. For Co > 0.62 the stresses 
at the equator are not saturated because the large scale 
flow is not fully developed. 
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3.3.2. Comparison to SOCA results 

To connect to earlier studies, the 3D simulation data and 
MTA-model results are compared to the analytical SOCA 
results of Kitchatinov & Riidiger ( 2005; see also Kitchatinov 
& Riidiger ll993p . Since analytical results are only available 
for the components relevant to the A-effect, only the Q xy 
and Q yz components are thus considered. Furthermore, we 
restrict the comparison to a subset of the models with Co = 
(0.15,1.5,5.4). 

The conventional way of writing the A-effect is (e.g. 
Riidiger HHEHI 

Q xy = u t Q, H cos 9 , (35) 



Q yz = v t D, Vsm6, (36) 

where z/ t = ^w rm s/fcf is the turbulent viscosity, and the 
dimensionless quantities H and V are given by 

H = H {1) sin 2 6 , (37) 
V = V {0) + V {1) cos 2 9 , (38) 

where V^°\ and — —H^ depend on the Coriolis 
number. 

The normalized stresses Qij / (vt^lo) are known from the 
simulations and the MTA-model, whereas H and V can be 
computed analytically for the turbulence model of KR05 
(see Appendix !b|. The results are shown in Fig.[l0]for three 
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Coriolis numbers. Our Coriolis number is smaller than that 
of KR05 by a factor of 2ir. 

For the horizontal stress, the SOCA results and the ro- 
tationally quenched MTA-model seem to fare similarly well. 
The former underestimates the magnitude for small Co and 
overestimates it for large Co, whereas for the latter the 
trend is exactly the opposite. If rotational isotropization is 
not taken into account, the agreement is poor for all Coriolis 
numbers considered here. 

For Q yz the standard MTA-model is almost spot on for 
Co = 0.15, but overestimates the magnitudes by at least a 
factor of two for the other cases. As discussed in the previ- 
ous section, the latitude distribution shows a mid-latitude 
maximum that is not present in the numerical simulations. 
When rotational isotropization is taken into account, at 
least the magnitude can be reconciled with the numerical 
results. The SOCA result does not fare very well in this 
case, predicting, in general, values that are too low and an 
incompatible latitude distribution with zero stress at the 
equator. 

3.4. Dependence on the amount of anisotropy 

Tabic [BT2l and Fig. [TT1 show the results for four sets of calcu- 
lations in which the Coriolis number is kept approximately 
constant at 0.3 whereas the turbulence anisotropy Ay is 
varied between —0.07 and —0.47. This is done by choosing 
suitable values for /o and /1 so that the rms- velocity stays 
approximately constant. 

From the figure it is seen that the ratio of the stress 
to the amount of anisotropy decreases monotonically as a 
function of Ay. Although part of the difference can be ex- 
plained by the somewhat smaller Co (0.31 as opposed to 
0.34 in the other cases) in the calculation with the largest 
Ay, the trend still persists. 

This trend can be understood as follows: from the ap- 
proximate relation Ay ~ 2r c Ay we obtain r c oc Ky/Ay. 
The decreasing trend of Q yz /Av oc Ay /Ay oc r c seen in 
the results suggests that the correlation time changes when 
the turbulence anisotropy is varied, i.e. r c decreases when 
Ay is increased. This is plausible since to change Ay differ- 
ent values of the forcing amplitudes, fo and /1 (see Eq.lTB")). 
need to be used resulting in differences in the turbulence. 

3.5. Dependence on Reynolds number 

Figure [T^] shows the off-diagonal stresses as functions of 
latitude and Reynolds number for a constant Co « 0.25 
(see also Table |R3)) . There is no a priori reason to expect 
that the stresses should depend on the Reynolds number if 
Re 3> 1 and if the turbulence anisotropy is kept constant. 
This is essentially what is borne out of the simulations, 
although there seems to be a weak decreasing trend as a 
function of Re for Re > 15 for Q xy and Q xz , although the 
results are within error bars. For Q yz , however, the largest 
Reynolds number case shows a distinct drop in comparison 
to the less turbulent cases. 

The decrease in the stresses as a function of the 
Reynolds number is likely to have the same origin as the 
decrease seen when the turbulence anisotropy is increased 
(see the previous section). In order to obtain the same Ay 
in the simulations different values, /o and /1 are needed 
for different values of v, the kinematic viscosity. More pre- 
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itude and vertical turbulence anisotropy Ay, see Eq. (j32|) . 
Line styles as indicated by the legend in the lower-most 
panel. Coriolis number in each case is roughly 0.3. 



cisely, the less the viscosity, the more difficult it is to obtain 
large anisotropy. Thus, for smaller v, higher ratio /1//0 is 
required to achieve a given Ay . Thus, if the interpretation 
that the higher the ratio /1//0, the smaller r c becomes is 
correct, the decreasing trend seen as a function of Re might 
be an artefact due to the small differences in the forcing be- 
tween the different runs. 
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3.6. Strouhal number from simulations 

The MTA-model reproduces the simulation results for the 
off-diagonal stresses reasonably well when St = 1 ... 2 is 
used. Now we turn to the numerical simulations to deter- 
mine St independently. Although we have not been able to 
provide any direct support of the basic MTA-assumption, 
Tij — —T~ 1 Qij 1 the MTA-model is still able to reproduce 
many of the features of the numerical simulations ade- 
quately. This implies that a value for St could be extracted 
from the MTA-relations. 

We consider two methods to determine the Strouhal 
number from the simulations: (i) MTA-relations derived for 



the Reynolds stresses, and (ii) similar relations for passive 
scalar transport under the influence of rotation. 

3.6.1. Off-diagonal Reynolds stresses versus MTA-relations 

Using the minimal tau-approximation and assuming a sta- 
tionary state where U = /force + /vise = 0, the off-diagonal 
Reynolds stresses can be derived from Eq. (|5J|, yielding 



Qxy — 2 QzT~xy(Qyy Qxx) 2 Qx^xy 
Qxz — 2 QzTxzQyz 2 QxTxzQxy j 
Qyz — 2 Q x T yz {Q zz Qyy 



2 Q z T yz Q xz : 



(39) 
(40) 
(41) 



where r,j are a set of relaxation times that allow for the pos- 
sibility that there can be different values of r for different 
components of the Reynolds stress. Here Qij are the com- 
ponents of the Reynolds stress tensor in the simulations. 
These equations can be solved for the via 

(42) 
(43) 



2 H 2 {Qyy Qxx) 2 QxQxz 
Qxz 

2 QzQyz 2 £2 X Q xy 
Qxy 

2 Q X (,Q zz Qyy) ~t~ ^ ^zQxz 



The first and the last of these equations, i.e. r xy and 
T yz, yield reasonable values for most cases, whereas t xz is 
less well-behaved, producing sign changes and showing no 
clear trend in magnitude as a function of latitude or rota- 
tion. This is because the two terms in the denominator of 
Eq. (|43|) tend to nearly cancel each other; see the upper- 
panels of Fig. [5] If, however, the values of r in Eqs. ([39]) 
to (|4"Tj) are considered the same, which is a basic MTA- 
assumption, it is possible to solve for Q xz in terms of the 
diagonal stresses 



4 Q x VI z t xz 

0' xz 



1 + 4 ttiT, 



{Qxx ~t" Q zz 2 Qyy ) , 



from which it follows that 



-1 1/2 



AQ x Q z {Qxx + Q zz ~ 2 Qyy) 4 Q^Qxz 



(45) 



(46) 



Although the assumption that all Tij are equal in Eqs. (f3"9"]) 
to (|41[) is not exactly realized in the numerical simulations, 
relation (I46[) gives results that are compatible with those 
from Eqs. (|42p and (|44p . The Strouhal number can now be 
computed from 



St 



(47) 



where are given by Eqs. (|42)) . (|46j) . and (|44)l . respectively. 
Representative results are given in Fig. 1131 

3.6.2. Passive scalar transport with rotation 

As an independent check of the dependence of the Strouhal 
number on rotation, we expand the passive scalar transport 
case, which was studied by Brandenburg et al. ( |2004p to 
cases where rotation is included. The numerical model is the 
same as in the runs presented so far, except that isotropic 
forcing is used with f = 0.01 - 0.03 and /i = in Eq. (II. 
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The turbulent passive scalar flux is denoted by T{ = UjC, 
where c is the fluctuation of passive scalar density, i.e. the 
passive scalar concentration per unit volume. Following the 
MTA-approach, we solve first for the time derivative 



*F i — ■ I ix^c . 



(48) 



where the fluctuation of the passive scalar field, neglecting 
diffusive terms, is given by 



- V • (uC + Uc + uc) 



(49) 



Here C and U are the mean passive scalar concentration 
and the mean velocity, respectively. Following Brandenburg 
et al. (|2004p we impose a large-scale gradient of passive 
scalar concentration according to VC = (0, 0, G) T . In what 
follows G = 0.1 is used. Now, assuming incompressibility 
and using the fact that U = 0, we arrive at 



.(2) 



T. 



(3) 



(50) 



where the last three terms denote the triple correlations. 
For the z-component of Eq. (|5"0)) the triple correlations are 
given by 

T z (1) = u z V ■ (uc) , Tj 2) = (uc) ■ Vu z , T z (3) = cVA (51) 
where h = c 2 In p is the reduced pressure (or enthalpy) . In 



the non-rotating case, + T^' — 0, and the contribu- 
tions from the momentum equation also cancel on average, 

(2) (3) 

i.e. Tz + T z — 0. The former relation follows from the 
periodic boundary conditions used and remains valid when 
rotation is added. The latter relation, however, is no longer 

(3) 

true and T z is now balanced by the Coriolis term. Thus, 



the MTA-assumption should be applied to T : 



(3) 



cV z h = t~ x T z 



(52) 



Assuming a stationary state in Eq. (|50p. the passive scalar 
fluxes can now be written as 



J y — TyQy Z G ~~ 2 £l Z TyJ~ X ~\~ 2 ^\ X Ty^F Z , 

J~ Z T Z Q ZZ G — 2 Q X T Z F y , 



(53) 
(54) 
(55) 



where we have retained the possibility that the values of t 
from different equations are unequal. 

In the passive scalar cases we consider isotropically 
forced turbulence for which Q xz ~ Q yz ~ even when 
rotation is included. Equations ([55]) (f55]) yield 



T x = ■ z = 



2 SX Z J~ y 



T 



2 £l x F z 2 Q Z J~ x 



(56) 
(57) 
(58) 



QzzG + 2 fl X J- y 

which can be used to compute the Strouhal numbers 

Stj = nu rms fcf, (59) 



where n are given by Eqs. (|56| to (|58|) . 

In the passive scalar case, the second and third order 
terms are indeed correlated (Brandenburg et al. I2004|) ac- 
cording to the basic MTA-assumption, and a Strouhal num- 
ber can be thus computed using 

St 3 = T 3 M rms fcf , (60) 

where 

T 3 =T Z /Wlh . (61) 

See Fig. [14] for representative results. 
3.6.3. Discussion 

For slow rotation, the Strouhal number is consistently be- 
tween one and three when computed from the Reynolds 
stresses (Fig. IT3)) . Similar values are obtained from the pas- 
sive scalar transport (Fig. I14[) when the Reynolds number 
is ten or larger. The Strouhal number computed from the 
triple correlations, St3, is more strongly dependent on the 
Reynolds number, but it seems to converge slowly towards 
a constant value near unity. These results are in line with 
the values required in the MTA-model and earlier studies in 
different contexts employing similar turbulence calculations 
(Brandenburg et al. 2004; Brandenburg & Subramanian 

I21KI51I2W1) . 

However, when the Coriolis number approaches or ex- 
ceeds unity, the Strouhal numbers computed from the equa- 
tions of the Reynolds stresses, Eqs. (|4"2")) - (|4*4")) . decrease 
rapidly so that for Co « 5 it has dropped at least by an or- 
der of magnitude (see Fig. [T3| . Similar results are obtained 
for the passive scalar transport under the influence of ro- 
tation, see Eqs. (f5"6")) - (|5"81) and the lower panels of Fig. Q3J 
The trend is clearer for St x and St y , whereas for St 2 the 
decreasing trend is seen only for rapid rotation, i.e. when 
Co > 2 ... 3. For slow rotation, however, St z is almost con- 
stant and increases somewhat when the Coriolis number 
approaches unity. These results seem to confirm the trend 
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seen earlier in convection simulations (Kapyla et al. 2005, 
I200fiap . 

The Strouhal number from the triple correlations fol- 
lows a trend similar to St z , with increasing values up to 
Co w 2 ... 3 after which there is a rapid decrease, see the 
uppermost panel of Fig. 1141 For low Reynolds number St 
can become negative in the range Co = 1 ... 5, hence the 
gaps in the corresponding data in Fig. 



4. Conclusions 

Turbulent momentum fluxes, which are described by the 
Reynolds stresses, were determined from numerical simula- 
tions of homogeneous rotating anisotropic turbulence. Since 
no large-scale shear is present, the generated Reynolds 
stresses correspond to contributions that are already 
present for uniform rotation. The resulting term is known as 
the A-effect (Krause & Riidiger ll974p . The component re- 
sponsible for the horizontal transport, Q xy , is positive and 
peaks around latitude 30° regardless of the Coriolis num- 
ber. The vertical component is predominantly negative and 
it always peaks at the equator. 

Although the numerical results for the A-effect broadly 
agree with analytical SOCA calculations (Kitchatinov & 
Riidiger 1993] I2005p . the MTA- model seems to reproduce 
certain features of the numerical results somewhat more 
closely. The present numerical results do not show the enig- 
matic results, such as the extreme latitude distribution of 
Q xy or a positive Q yz for rapid rotation, which have been re- 
ported from convection simulations (e.g. Chan 2001; Kapyla 
et al. 2004). The difference lies most likely in our neglecting 
stratification and heat fluxes. The exact manner in which 
they affect the Reynolds stresses is not within the scope of 
the present paper, but should be investigated more closely 
in the future. 

By applying the minimal tau approximation closure re- 
lation to the Reynolds stress equation, qualitatively similar 
results are obtained, but the magnitude of the stresses is in 
general too large. The vertical flux in the MTA-model, how- 
ever, has a maximum at mid-latitudes for intermediate and 
rapid rotation. Adding an empirical rotational isotropiza- 
tion term (motivated in Sect. 13. lj) also brings the magnitude 
in line with the 3D simulations. Although adding this term 
with this particular form has no rigorous theoretical basis, 
we can see that phenomenological effects of isotropization 
of turbulence due to rotation are indeed seen in the simu- 
lations and that the term is thus justified. 

Another drawback of the MTA-model is that the diago- 
nal components of the Reynolds tensor are rather badly 
reproduced since the nonlinear effects of rotation mani- 
fest in the numerical simulations are not explicitly taken 
into account. The empirically added rotational isotropiza- 
tion term augments the magnitudes, but not the latitude 
distribution. Furthermore, no direct evidence of the validity 
of the MTA-assumption Qij = —Tij/r was found in the nu- 
merical simulations. Contrasting the behavior of the diago- 
nal components to the fairly good correspondence between 
the numerical simulations and the MTA-model for the off- 
diagonal components leads us to conclude that, where the 
behavior of the diagonal components is dominated by the 
inadequately modeled nonlinear effects, the off-diagonals 
are fairly well presented by the linear terms. 

A Strouhal number of order unity in the MTA-model 
gives best fits to the numerical results. Fitting the numer- 
ical results to expressions derived under the MTA, similar 
values of St are found for slow rotation. For Coriolis num- 
bers of order unity or larger, however, the Strouhal number 
obtained in this manner decreases rapidly. In the passive 
scalar case, the situation is somewhat more complex, al- 
though a similar decreasing trend of the Strouhal number 
is recovered for rapid rotation, see Fig. [14] These results are 
in accordance with earlier results from convection Simula- 
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tions (Kapyla et al. 120051 12006a) using Reynolds stresses or 
correlation analysis of the velocity field. 

A related aspect in turbulent transport that requires 
closer study is the turbulent viscosity (see preliminary re- 
sults in Kapyla & Brandenburg 12007)) and the possibility 
of a A-effect due to the anisotropy induced by a large-scale 
shear flow (Leprovost & Kim 2007) . These matters will be 
considered in more detail in a future publication. 
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Appendix A: Stationary solutions for the stresses 
from MTA 

Setting Qij = and using the MTA-closure TJy = —Qij/r, 
and parameterizing the contributions of the forcing as 
Uifj + Ujfi = Qj?V T m Eq. (|2"2"]) gives six equations for 
the six unknowns of the symmetric tensor Qij, i.e. 

<°) (A.l) 

(A.2) 
(A.3) 

Q yy = 4 Q x tQ vz - 4 n z rQ xy + , (A.4) 

Qyz — 2 Q x T~(yQ zz Qyy) ^ QzTQxz ; (A. 5) 

Qzz = -4Q x TQ yz + Q<£} . (A.6) 
It is possible to solve for Qij in terms of r, Q x , Q z , and 
Qlf. From Eqs. (|A.1|) to (|A.6[) it is clear that only three 
components of Qij are independent. Thus it is sufficient to 
solve for the off-diagonal components. After some algebra 
we arrive at 



Qxx — 4 O z7~Q xy + Q X x ' 
Qxy — ^ Qz7~{Qyy Qxx) + ^ £I x tQ q 
Qxz 2 Qz^~Qyz ^ £1 X 7~Q xy j 

)(0) 



Qxy = 2n,r(j| , )«&>-Q«S) 



24^o zT 3 Xo(Q (o)_ Q (o) )! (A7) 

Qxz — 2 Q>z7~Qyz 2 Q>x7~Q xy i (-^■■^) 
Qyz = 20 x tAT 1 (1 + 144r»^2 r 4 Ko)(Q (0) _ g(0)j + 

24n x nlr 3 K (Q$-QW), (A.9) 

where 

1 

K n = 



1 + 20^7-2 +64 0^4 



1 + 5Co MTA + 8Co MTA 



(A.10) 



1 +4 02 T 2 + 16 02 r 2 



1 



1 + Comta(1 + 3 sin 2 i 
and Comta is given by Eq. (|2"7) . 

In the present study the forcing is such that Q 



,(A.ll) 



(0) 
yy 



Qxx = 0, so the equations reduce to 

Qxy = 24O2O J2 r 3 A (gi°;-Q^), 

Qxz — 2 QzfQyz 2 Qx7~Qxy ) 
Qyz 



(A.12) 
(A.13) 

2fi x T# 1 (l + 144O^r 4 A )(Qi z ) - Q(°A (A.14) 



Appendix B: Coefficients of the A-effect from 
SOCA 

The fluxes of angular momentum have commonly been pa- 
rameterized by Eqs. (|35p and (|36[) . and the normalized 
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fluxes by Eqs. d2ZJ) and ||35J|. Kitchatinov & Riidiger (j2"UU5|) 
computed these coefficients using SOCA 

H = ff (1) (^)sin 2 6l , (B.l) 
V = (0*) - (fi*) cos 2 6> , (B.2) 

where fi* = 2 f2o T to is their definition of the Coriolis num- 
ber, and Tto the turnover time. Note that there is a differ- 
ence of 2-7T in comparison to our definition, Eq. fll7[) . For 
simplicity, we retain f2* in the expressions that follow. The 
coefficients H^ 1 ' and are given by 

Hil) = [Ji(n*) +0^(0*)] , (B.3) 

^ = (ip) [ J o(^)+a/oM, (B.4) 

where i corr is the correlation length, i? p the density scale 
height, and a = 2 an 'anisotropy parameter' that reduces 
the amount of anisotropy for slow rotation. 
The functions Jj and are given by 

7 » ^(- 19 -T^ + i T" cto ' in «)' (B ' 5) 
7 ' = ^(- 15 -T^ + T» cto ''4 (&6) 

1 ( 2 2 ^2 + g \ 



2ftfV" ' ""* 

- 12 ft? -45 



arctanft* I . (B 



In our model there is no density stratification, so l CO rr/H p 
is taken to be a free parameter that we choose to be equal 
to unity. 
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Table B.l. Summary of the turbulence anisotropies and normalized Reynolds stresses, Qij = Qij/u^ ms , for the Re 
12 ... 14 calculations, fcf/fci = 5, v = 2 ■ 10~ 3 , and grid resolution 64 3 was used in all runs. 



Run 




6 


Co 


^rms 


Qxx 


Qyy 


Qzz 




A v 


Qxy 


Qxz 


Qyz 


64al 


0.02 


0° 


0.06 


0.127 


0.236 


0.240 


0.525 


0.002 


-0.571 


-0.001 


-0.000 


0.003 


64a2 


0.02 


15° 


0.06 


0.127 


0.239 


0.242 


0.523 


0.003 


-0.567 


-0.001 


-0.002 


-0.004 


64a3 


0.02 


30° 


0.06 


0.127 


0.237 


0.243 


0.523 


0.004 


-0.564 


-0.001 


-0.003 


-0.001 


64a4 


0.02 


45° 


0.06 


0.127 


0.238 


0.244 


0.522 


0.005 


-0.561 


-0.000 


-0.004 


-0.014 


64a5 


0.02 


60° 


0.06 


0.127 


0.237 


0.246 


0.521 


0.008 


-0.559 


-0.000 


-0.003 


-0.018 


64a6 


0.02 


75° 


0.06 


0.127 


0.237 


0.246 


0.521 


0.010 


-0.559 


-0.001 


-0.002 


-0.021 


64a7 


0.02 


90° 


0.06 


0.127 


0.236 


0.247 


0.521 


0.010 


-0.558 


-0.001 


-0.001 


-0.021 



64a8 


0.05 


0° 


0.15 


0.127 


0.244 


0.246 


0.514 


0.003 


-0.537 


0.000 


0.001 


0.002 


64a9 


0.05 


15° 


0.15 


0.127 


0.244 


0.248 


0.512 


0.003 


-0.532 


-0.000 


-0.005 


-0.007 


64al0 


0.05 


30° 


0.15 


0.127 


0.243 


0.250 


0.511 


0.006 


-0.529 


0.002 


-0.010 


-0.018 


64all 


0.05 


45° 


0.15 


0.127 


0.242 


0.253 


0.509 


0.011 


-0.522 


0.004 


-0.012 


-0.027 


64al2 


0.05 


60° 


0.15 


0.127 


0.239 


0.258 


0.507 


0.019 


-0.516 


0.005 


-0.009 


-0.035 


64al3 


0.05 


75° 


0.15 


0.127 


0.236 


0.263 


0.505 


0.027 


-0.512 


0.004 


-0.005 


-0.040 


64al4 


0.05 


90° 


0.15 


0.127 


0.234 


0.265 


0.505 


0.032 


-0.510 


-0.000 


0.000 


-0.042 



64al5 


0.1 


0° 


0.31 


0.126 


0.255 


0.258 


0.490 


0.003 


-0.468 


-0.000 


0.002 


0.002 


64al6 


0.1 


15° 


0.31 


0.127 


0.256 


0.257 


0.490 


0.001 


-0.468 


0.001 


-0.009 


-0.004 


64al7 


0.1 


30° 


0.31 


0.127 


0.257 


0.257 


0.490 


-0.000 


-0.466 


0.005 


-0.019 


-0.014 


64al8 


0.1 


45° 


0.31 


0.127 


0.254 


0.261 


0.489 


0.007 


-0.462 


0.011 


-0.022 


-0.028 


64al9 


0.1 


60° 


0.31 


0.127 


0.247 


0.272 


0.484 


0.025 


-0.448 


0.014 


-0.017 


-0.039 


64a20 


0.1 


75° 


0.31 


0.127 


0.239 


0.284 


0.480 


0.045 


-0.437 


0.010 


-0.008 


-0.045 


64a21 


0.1 


90° 


0.31 


0.127 


0.235 


0.289 


0.479 


0.054 


-0.434 


-0.000 


0.000 


-0.047 



64a22 


0.2 


0° 


0.62 


0.127 


0.278 


0.280 


0.446 


0.002 


-0.333 


-0.001 


0.001 


-0.000 


64a23 


0.2 


15° 


0.62 


0.127 


0.277 


0.276 


0.451 


-0.001 


-0.350 


0.001 


-0.008 


0.005 


64a24 


0.2 


30° 


0.62 


0.127 


0.276 


0.268 


0.460 


-0.008 


-0.376 


0.005 


-0.021 


-0.002 


64a25 


0.2 


45° 


0.62 


0.127 


0.273 


0.268 


0.462 


-0.005 


-0.382 


0.013 


-0.026 


-0.015 


64a26 


0.2 


60° 


0.61 


0.128 


0.267 


0.280 


0.457 


0.012 


-0.367 


0.019 


-0.021 


-0.027 


64a27 


0.2 


75° 


0.61 


0.128 


0.255 


0.298 


0.451 


0.043 


-0.350 


0.013 


-0.010 


-0.035 


64a28 


0.2 


90° 


0.61 


0.128 


0.251 


0.303 


0.450 


0.053 


-0.345 


-0.001 


-0.000 


-0.036 



64a29 


0.5 


0° 


1.47 


0.133 


0.303 


0.301 


0.400 


-0.003 


-0.197 


-0.001 


0.003 


0.000 


64a30 


0.5 


15° 


1.48 


0.132 


0.299 


0.295 


0.410 


-0.004 


-0.225 


0.000 


-0.009 


0.005 


64a31 


0.5 


30° 


1.49 


0.131 


0.296 


0.286 


0.422 


-0.010 


-0.262 


0.003 


-0.019 


0.005 


64a32 


0.5 


45° 


1.50 


0.131 


0.290 


0.283 


0.431 


-0.007 


-0.288 


0.008 


-0.027 


-0.004 


64a33 


0.5 


60° 


1.51 


0.130 


0.286 


0.289 


0.430 


0.003 


-0.285 


0.012 


-0.022 


-0.013 


64a34 


0.5 


75° 


1.50 


0.130 


0.279 


0.303 


0.422 


0.024 


-0.262 


0.007 


-0.011 


-0.018 


64a35 


0.5 


90° 


1.50 


0.131 


0.274 


0.305 


0.422 


0.035 


-0.261 


-0.001 


-0.000 


-0.018 



64a36 


1.0 


0° 


2.78 


0.141 


0.310 


0.310 


0.384 


-0.000 


-0.148 


0.001 


0.003 


0.001 


64a37 


1.0 


15° 


2.80 


0.140 


0.306 


0.308 


0.390 


0.002 


-0.166 


0.000 


-0.008 


0.003 


64a38 


1.0 


30° 


2.84 


0.138 


0.303 


0.296 


0.405 


-0.007 


-0.212 


0.001 


-0.004 


0.002 


64a39 


1.0 


45° 


2.88 


0.136 


0.292 


0.299 


0.412 


0.007 


-0.233 


0.004 


-0.019 


-0.003 


64a40 


1.0 


60° 


2.90 


0.135 


0.288 


0.302 


0.414 


0.014 


-0.238 


0.007 


-0.012 


-0.006 


64a41 


1.0 


75° 


2.89 


0.136 


0.274 


0.316 


0.414 


0.042 


-0.239 


0.006 


-0.004 


-0.009 


64a42 


1.0 


90° 


2.87 


0.137 


0.256 


0.320 


0.428 


0.065 


-0.279 


0.000 


0.002 


-0.009 



64a43 


2.0 


0° 


5.26 


0.149 


0.293 


0.294 


0.417 


0.001 


-0.247 


-0.001 


0.001 


0.001 


64a44 


2.0 


15° 


5.27 


0.149 


0.298 


0.300 


0.406 


0.002 


-0.214 


0.001 


-0.028 


-0.001 


64a45 


2.0 


30° 


5.41 


0.145 


0.311 


0.310 


0.384 


-0.001 


-0.147 


0.002 


-0.030 


-0.003 


64a46 


2.0 


45° 


5.41 


0.144 


0.293 


0.318 


0.393 


0.024 


-0.174 


0.003 


-0.014 


-0.004 


64a47 


2.0 


60° 


5.46 


0.144 


0.264 


0.337 


0.404 


0.073 


-0.206 


0.005 


0.010 


-0.006 


64a48 


2.0 


75° 


5.46 


0.144 


0.228 


0.350 


0.427 


0.122 


-0.277 


0.006 


0.012 


-0.005 


64a48.5 


2.0 


82.5° 


5.44 


0.144 


0.220 


0.348 


0.437 


0.128 


-0.306 


0.002 


0.007 


-0.007 


64a49 


2.0 


90° 


5.43 


0.145 


0.200 


0.361 


0.444 


0.161 


-0.329 


0.000 


-0.000 


-0.005 
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Table B.2. Summary of the turbulence anisotropies and normalized Reynolds stresses, = Qij/u^ ms , for a set of runs 
with varying turbulence anisotropy. Oq = 0.1, fcf/fci = 5, v — 2 • 10~ 3 , and grid resolution 64 3 was used in all runs. 



Run 


n 


e 


Co 


li rm s 


Qxx 


Qyy 


Qzx 


An 


Ay 


Qxy 


Qxz 


Qyz 


64al5 


0.1 


0° 


0.31 


0.126 


0.255 


0.258 


0.490 


0.003 


-0.468 


-0.000 


0.002 


0.002 


64al6 


0.1 


15° 


0.31 


0.127 


0.256 


0.257 


0.490 


0.001 


-0.468 


0.001 


-0.009 


-0.004 


64al7 


0.1 


30° 


0.31 


0.127 


0.257 


0.257 


0.490 


-0.000 


-0.466 


0.005 


-0.019 


-0.014 


64al8 


0.1 


45° 


0.31 


0.127 


0.254 


0.261 


0.489 


0.007 


-0.462 


0.011 


-0.022 


-0.028 


64al9 


0.1 


60° 


0.31 


0.127 


0.247 


0.272 


0.484 


0.025 


-0.448 


0.014 


-0.017 


-0.039 


64a20 


0.1 


75° 


0.31 


0.127 


0.239 


0.284 


0.480 


0.045 


-0.437 


0.010 


-0.008 


-0.045 


64a21 


0.1 


90° 


0.31 


0.127 


0.235 


0.289 


0.479 


0.054 


-0.434 


-0.000 


0.000 


-0.047 



64cl5 


0.1 


0° 


0.34 


0.115 


0.275 


0.275 


0.453 


0.000 


-0.357 


0.001 


0.003 


0.000 


64cl6 


0.1 


15° 


0.34 


0.116 


0.275 


0.276 


0.454 


0.001 


-0.357 


0.003 


-0.009 


-0.007 


64cl7 


0.1 


30° 


0.34 


0.116 


0.277 


0.275 


0.452 


-0.002 


-0.353 


0.006 


-0.017 


-0.018 


64cl8 


0.1 


45° 


0.34 


0.116 


0.275 


0.280 


0.449 


0.006 


-0.343 


0.011 


-0.019 


-0.030 


64cl9 


0.1 


60° 


0.34 


0.115 


0.271 


0.290 


0.443 


0.019 


-0.324 


0.014 


-0.015 


-0.038 


64c20 


0.1 


75° 


0.34 


0.116 


0.265 


0.300 


0.439 


0.035 


-0.314 


0.010 


-0.007 


-0.043 


64c21 


0.1 


90° 


0.34 


0.116 


0.262 


0.304 


0.439 


0.042 


-0.311 


0.002 


0.001 


-0.044 



64b 15 


0.1 


0° 


0.33 


0.119 


0.299 


0.299 


0.406 


0.000 


-0.215 


0.001 


0.001 


0.002 


64b 16 


0.1 


15° 


0.33 


0.119 


0.300 


0.299 


0.405 


-0.000 


-0.211 


0.002 


-0.005 


-0.005 


64bl7 


0.1 


30° 


0.33 


0.119 


0.302 


0.299 


0.404 


-0.003 


-0.207 


0.004 


-0.011 


-0.012 


64bl8 


0.1 


45° 


0.33 


0.119 


0.301 


0.303 


0.400 


0.002 


-0.195 


0.007 


-0.012 


-0.019 


64b 19 


0.1 


60° 


0.33 


0.119 


0.300 


0.308 


0.396 


0.008 


-0.185 


0.010 


-0.011 


-0.024 


64b20 


0.1 


75° 


0.33 


0.119 


0.297 


0.313 


0.395 


0.016 


-0.179 


0.006 


-0.004 


-0.026 


64b21 


0.1 


90° 


0.33 


0.119 


0.295 


0.317 


0.393 


0.022 


-0.175 


0.003 


-0.001 


-0.027 



64dl5 


0.1 


0° 


0.34 


0.115 


0.321 


0.318 


0.365 


-0.002 


-0.091 


0.001 


0.000 


0.001 


64dl6 


0.1 


15° 


0.34 


0.115 


0.320 


0.320 


0.364 


0.000 


-0.087 


0.001 


-0.004 


-0.002 


64dl7 


0.1 


30° 


0.34 


0.115 


0.321 


0.320 


0.363 


-0.001 


-0.086 


0.001 


-0.006 


-0.005 


64dl8 


0.1 


45° 


0.34 


0.115 


0.321 


0.323 


0.360 


0.002 


-0.076 


0.003 


-0.006 


-0.009 


64dl9 


0.1 


60° 


0.34 


0.115 


0.324 


0.323 


0.357 


-0.001 


-0.066 


0.004 


-0.006 


-0.012 


64d20 


0.1 


75° 


0.34 


0.115 


0.323 


0.325 


0.356 


0.002 


-0.065 


0.003 


-0.004 


-0.013 


64d21 


0.1 


90° 


0.34 


0.115 


0.319 


0.327 


0.357 


0.009 


-0.069 


0.000 


-0.000 


-0.013 
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Table B.3. Turbulence anisotropies and Reynolds stresses for varying Reynolds numbers, fcf/fci = 5 and £1 = 0.1 were 
used in all runs. Qij = Qij/u^ 



,2 

rms ' 



Run 


Grid 


e 


Co 


!/[10~ 3 


] Re 


Urms 




Qyy 


Qzz 


An 


Av 


Qxy 


Qxz 


Qyz 


32cl5 


32 3 


0° 


0.255 


5 


6.0 


0.154 


0.274 


0.273 


0.454 


-0.000 


-0.361 


0.001 


-0.001 


-0.001 


32cl6 


32 3 


15° 


0.255 


5 


6.0 


0.154 


0.274 


0.273 


0.454 


-0.001 


-0.361 


0.002 


-0.007 


-0.009 


32cl7 


32 3 


30° 


0.255 


5 


6.0 


0.154 


0.273 


0.275 


0.453 


0.001 


-0.359 


0.004 


-0.011 


-0.018 


32cl8 


32 3 


45° 


0.255 


5 


6.0 


0.154 


0.272 


0.278 


0.451 


0.006 


-0.352 


0.006 


-0.012 


-0.027 


32cl9 


32 3 


60° 


0.255 


5 


6.0 


0.154 


0.270 


0.283 


0.448 


0.014 


-0.344 


0.007 


-0.010 


-0.034 


32c20 


32 3 


75° 


0.255 


5 


6.0 


0.154 


0.267 


0.288 


0.446 


0.021 


-0.338 


0.006 


-0.006 


-0.039 


32c21 


32 3 


90° 


0.255 


5 


6.0 


0.154 


0.266 


0.290 


0.445 


0.025 


-0.335 


0.002 


-0.001 


-0.041 


64el5 


64 3 


0° 


0.254 


2 


15 


0.155 


0.273 


0.274 


0.454 


0.001 


-0.361 


0.000 


0.000 


0.001 


64el6 


64 3 


15° 


0.254 


2 


15 


0.155 


0.273 


0.274 


0.454 


0.001 


-0.361 


0.000 


-0.009 


-0.008 


64el7 


64 3 


30° 


0.254 


2 


15 


0.155 


0.273 


0.276 


0.451 


0.003 


-0.353 


0.004 


-0.016 


-0.019 


64el8 


64 3 


45° 


0.254 


2 


15 


0.155 


0.271 


0.283 


0.446 


0.011 


-0.339 


0.009 


-0.017 


-0.030 


64el9 


64 3 


60° 


0.254 


2 


15 


0.155 


0.269 


0.290 


0.442 


0.021 


-0.324 


0.010 


-0.013 


-0.036 


64e20 


64 3 


75° 


0.253 


2 


15 


0.155 


0.265 


0.297 


0.439 


0.032 


-0.316 


0.006 


-0.007 


-0.041 


64e21 


64 3 


90° 


0.253 


2 


15 


0.155 


0.263 


0.300 


0.437 


0.037 


-0.311 


-0.001 


0.000 


-0.041 



128bl5 


128 3 


0° 


0.252 


1 


31 


0.156 


0.274 


0.275 


0.452 


0.000 


-0.355 


0.001 


-0.001 


0.002 


128bl6 


128 3 


15° 


0.251 


1 


31 


0.156 


0.276 


0.275 


0.450 


-0.001 


-0.349 


0.001 


-0.011 


-0.007 


128bl7 


128 3 


30° 


0.251 


1 


31 


0.156 


0.278 


0.277 


0.447 


-0.001 


-0.339 


0.004 


-0.016 


-0.017 


128bl8 


128 3 


45° 


0.251 


1 


31 


0.156 


0.278 


0.283 


0.440 


0.006 


-0.320 


0.009 


-0.019 


-0.026 


128bl9 


128 3 


60° 


0.250 


1 


31 


0.157 


0.275 


0.291 


0.435 


0.016 


-0.304 


0.011 


-0.015 


-0.035 


128b20 


128 3 


75° 


0.250 


1 


31 


0.157 


0.270 


0.299 


0.432 


0.029 


-0.295 


0.007 


-0.007 


-0.037 


128b21 


128 3 


90° 


0.250 


1 


31 


0.157 


0.268 


0.303 


0.430 


0.035 


-0.289 


-0.000 


-0.001 


-0.038 


256al5 


256 3 


0° 


0.251 


0.5 


61 


0.156 


0.275 


0.276 


0.449 


0.001 


-0.347 


-0.000 


-0.000 


0.003 


256al6 


256 3 


15° 


0.252 


0.5 


61 


0.156 


0.277 


0.276 


0.447 


-0.002 


-0.341 


0.001 


-0.008 


-0.005 


256al7 


256 3 


30° 


0.252 


0.5 


61 


0.156 


0.276 


0.277 


0.448 


0.001 


-0.342 


0.004 


-0.015 


-0.011 


256al8 


256 3 


45° 


0.250 


0.5 


61 


0.157 


0.275 


0.282 


0.443 


0.006 


-0.330 


0.009 


-0.015 


-0.020 


256al9 


256 3 


60° 


0.251 


0.5 


61 


0.156 


0.273 


0.289 


0.438 


0.016 


-0.314 


0.009 


-0.013 


-0.027 


256a20 


256 3 


75° 


0.250 


0.5 


61 


0.157 


0.269 


0.296 


0.435 


0.027 


-0.306 


0.006 


-0.007 


-0.030 


256a21 


256 3 


90° 


0.250 


0.5 


61 


0.157 


0.264 


0.301 


0.435 


0.037 


-0.305 


0.000 


-0.002 


-0.032 



